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1.  Introduction 

The  low  chromosphere’s  heating  by  the  series  of  shock  waves  is  investigated.  Waves 
are  initially  generated  in  the  solar  photosphere.  Due  to  propagation  in  the  strong  stratified  solar 
atmosphere  the  waves  are  overturning  to  the  shock  waves.  The  goal  of  present  paper  is  to  inves¬ 
tigate  it  using  2D  approximation. 

To  describe  the  solar  atmosphere  one-fluid  hydrodynamics  with  taken  into  consideration 
radiative  heat  transfer  and  absence  of  the  thermodynamic  equilibrium  is  used.  Computer  model¬ 
ing  is  performed  for  one-  and  two-dimensional  formulation  of  problem.  The  results  obtained  are 
important  for  understanding  the  physics  of  star  atmospheres’  reaction  on  the  dynamical  distur¬ 
bances  like  flares,  mass  surges  and  heating  by  penetrating  convection. 

2.  Model  approximation 

Solar  atmosphere  is  mixed  hydrogen-helium  plasma  in  strong  gravitational  field  (for  exact 
chemical  compositions  see  [1]).  The  hydrodynamic  approximation  is  proved  to  be  valid  due  to 
smallness  of  particle’s  free  path  compared  with  characteristic  spatial  scale  of  the  atmosphere 

(scale  of  height  A  =  p/(pg)  =  {d\n(p)j  d  r\l )  and  smallness  oftypical  time  of  the  kinetic  pro¬ 
cesses  (like  ionization,  recombination,  collision  rate  and  so  on)  compared  with  time  scale  of  the 
investigated  processes.  Thermoconductivity  should  be  taken  into  consideration  because  thermo¬ 
conductivity  coefficient  depends  only  on  temperature  and  thermoconductivity  becomes  impor¬ 
tant  in  the  chromosphere’s  upper  layer,  where  energy’s  volume  density  becomes  small  due  to  a 
strong  stratification.  Instead  of  thermoconductivity,  viscosity  is  negligible  in  all  bulk  of  photo¬ 
sphere  and  chromosphere. 

The  radiative  transfer  is  the  most  important  source  of  energy  loss  in  solar  atmosphere. 
Photon’s  free  path  even  in  the  middle  chromosphere  significantly  exceeds  the  height’s  scale,  so 
typical  value  of  optical  density  is  infinitesimal.  Optical  density  becomes  comparable  with  unity 
only  in  the  photosphere  and  it  proves  the  using  of  volumetric  loss  approximation  for  radiative 
transfer  describing  [2,  3  and4]. 

Solar  atmosphere  is  not  in  thermodynamic  equilibrium  (TE),  so  ions’  and  electrons’ 
densities,  absorption  factor  and  emittance  are  founded  by  solving  the  system  of  kinetic  equa¬ 
tions  included  the  all  significant  kinds  of  ionisation  and  recombination  processes.  The  main 
ionization  process  is  ionization  by  the  electron-atom  collision,  and  main  recombination  process 
is  free-bound  recombination  instead  of  reversal  process  of  the  tree  body  recombination.  That 
kind  of  equilibrium  is  typical  for  stellar  atmospheres  and  named  as  corona  equilibrium.  The 
densities  of  ions,  electrons  and  neutral  atoms  for  solar  chromosphere  are  founded  in  [2,  3]  as 
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Figure  1.  Dependence  of  ionization  parameter  on  Figure  2.  Dependence  of  heat  conductivity  coeffi- 
temperature.  cient  on  temperature 


well,  as  volumetric  loss  power  per  volume  unit.  (Rosseland  opacities  (absorption  factor,  aver¬ 
aged  over  spectra)  usually  are  founded  in  case  of  thermodynamic  equilibrium  [5,  references 
herein].  They  are  useless  because  of  the  correlation  between  absorption  factor  and  emittance  is 
consequence  of  thermodynamic  equilibrium  [6]). 

So  a  two-dimensional  system  of  equations  which  describes  the  solar  atmosphere’s  plasma 
in  hydrodynamic  approximation  with  radiative  heat  transfer  and  thermoconductivity  considered  is: 
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Here  p  designates  plasmas  density,  p  and  s  =  p/{{y  -1)  •  p)  —  denote  pressure  and  mass  den¬ 
sity  of  plasma  internal  energy,  w  =  ps  +  p(u2  +  v2)/ 2,  u,v  -  components  of  speed, 
g  =  274  m/s2  -  acceleration  of  gravity  (x  axis  is  antiparallel  to  the  g  vector),  %  =  %(T)—  coef¬ 
ficient  of  thermoconductivity,  Prad  ( p.  T)  -  volumetric  loss  power  caused  by  emittance 
(founded  in  [2-4]:  Prad(p,T)  =  nenH  ■ Lr(T ) ).  In  case  of  corona  equilibrium  ne  =  a(T)  ■  nH  , 
so  volumetric  loss  power  is  proportional  to  the  plasma  density  raised  to  the  second  power. 
Functions  %(T),a(T),Prad  (p,T)/(nenH )  are  plotted  on  Figs.  1-3. 

The  system  of  equations  is  solved  in  rectangular  area.  On  an  axis  "y"  the  periodic  bound¬ 
ary  conditions  are  set.  On  a  low  bound  the  speed  sets  as  a  function  from  time  and  coordinate 
(simulating  a  source  of  acoustic  waves).  On  a  upper  bound  the  boundary  condition  of  waves’ 
issue  is  putted.  All  equations  are  dimensionless  and  the  units  of  physical  parameters  are  listed 
below: 

[r]  =  107  cm  [?]  =  101’5  s  [w]  =  [r]/[t]  [p\  =  1  (T6  gr/cm3 

[Z]  =  104K  \p]  =  105  erg/cm3  [s]  =  1011  erg/gr  [g]  = 
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on  temperature. 


eral  times  for  one-dimensional  heating  (units  are 


dimensionless). 


Figure  6.  Temperature  distributions  over  height 


several  times  for  one-dimensional  heating  (units 
are  dimensionless) 


for  several  times  for  one-dimensional  heating 
(units  are  dimensionless) 


3.  Atmosphere’s  heating  by  the  series  of  acoustic  waves 

Earlier  the  atmosphere’s  heating  was  investigated  in  one-dimensional  formulations  of  a  prob¬ 
lem  and  numerical  scheme’s  realization  together  with  model’s  approximations  force  the  height’s 
range  to  be  less  than  2500  km  [7,  8].  The  solar  atmosphere  is  hardly  inhomogeneous  in  all  directions, 
so  the  research  of  a  problem  at  least  in  two-dimensional  formulation  is  of  strong  interest. 

With  the  purpose  to  find  out  influence  of  spatial  inhomogeneity  on  atmosphere’s  heating 
two  formulations  of  problem  are  performed.  In  the  first  computer  simulation  the  wave’s  source 
does  not  depend  on  “j/”  coordinate,  so  the  formulation  is  one-dimensional.  On  the  left  boundary 
speed  is  set  to  u(y,t)  =  Uq  sin(cof) ,  where  u0  =100  ml s  ,  2n / a>  =  100  s  are  the  typical  pa¬ 
rameters  of  acoustic  waves  in  the  solar  photosphere  [9].  Weak  acoustic  waves,  initially  gener¬ 
ated  on  the  low  boundary,  are  propagating  to  the  upper  boundary,  so  waves’  amplitude  is  expo¬ 
nentially  increasing  with  height  due  to  plasmas’  stratification.  If  speed  of  flow  becomes  bigger 


5 


than  sound  speed,  waves  are  turning  into  shock  waves’  series  (Fig.  4).  The  result  of  heating  by 
shock  waves  is  a  significant  decreasing  of  height’s  scale  (Fig.  5),  so  after  5-8  shock  wave’s 
transitions  the  atmosphere  becomes  quasi-stationary  and  heating  by  shock  waves  is  compen¬ 
sated  by  radiation  losses.  When  heating  process  relaxes  to  stationary  state,  the  shock  wave  am¬ 
plitude  becomes  of  order  1.26- 106  cm/s.  Distribution  of  plasmas’  parameters  (for  corresponding 
parameters  of  wave’s  source)  are  showed  in  Figs.  4,  5,  and  6. 

In  the  second  computer  simulation  wave’s  source  is  inhomogeneous  and  his  spatial  scale 
is  A  =  5.0  -108  cm  :  u(y,t)  =  Uq  •  sin(cof)  -sm(2ny/X) .  The  initial  atmosphere  was  taken  to  be 
heated  by  homogeneous  wave’s  source  instead  of  be  in  isothermal  hydrostatic  equilibrium. 
Though  the  wave’s  source’s  scale  is  more  than  100  times  bigger  than  scale  of  height,  influence 
of  that  inhomogeneity  is  sufficient.  Shock  waves  get  curved  and  their  fronts’  expansion  reduces 
the  increasing  of  wave’s  amplitude  with  height  as  well,  as  power  of  atmosphere  heating.  In  do¬ 
mains  upper  than  15000  km  shock  wave  develops  into  a  strong  curved  arc,  which  vertical  size  is 
of  order  of  scale  height  and  horizontal  size  is  2-4  times  greater  than  X.  It  leads  to  elimination  of 
waves’  series  regular  structure  and  explains  the  absence  of  well-defined  shock  waves’  series  in 
the  solar  corona  according  to  observational  data.  Minimal  and  maximal  values  of  the  some 
plasma  flow  parameters  during  period  of  wave’s  source  is  showed  on  figures  7-10.  On  figures 
11,  12  typical  distributions  of  the  temperature  and  vertical  component  of  speed  are  showed. 
Parameters  u0 ,  co  are  the  same  as  for  previous  one-dimensional  variant. 

4.  Conclusion 

In  present  paper  the  solar  atmosphere’s  heating  is  investigated  using  one-dimensional  and 
two-dimensional  formulation  of  the  problem.  It  is  proved  that  in  strong  stratified  medium  the 
influence  of  wave’s  source’s  inhomogeneity  becomes  important  on  heights,  which  are  compa¬ 
rable  with  several  scales  of  height  instead  of  scale  of  source’s  inhomogeneity.  Performed  inves¬ 
tigation  limits  the  earlier  obtained  one-dimensional  result’s  applicability  to  the  domains  with 
heights  of  order  700  -  1200  km,  because  of  typical  wave  source  scale  is  700-1500  km  (granule) 
or  3000  -7000  km  (length  of  coherence  of  5-min  oscillations). 
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Figure  7.  Minimal  over  period  temperature  for 
atmosphere,  heated  by  non-uniform  wave’s 
source  (units  are  dimensionless). 
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Figure  8.  Maximal  over  period  temperature  for 
atmosphere,  heated  by  non-unifonn  wave’s 
source  (units  are  dimensionless). 
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Figure  9.  Minimal  over  period  vertical  speed’s 
component  for  atmosphere,  heated  by  non-uniform 


wave’s  source  (units  are  dimensionless) 


Figure  1 1 .  Spatial  distribution  of  the  vertical  speed’s 
component  for  the  atmosphere,  heated  by  non- 
uniform  waves’  source  (all  units  are  dimensionless) 
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Figure  10.  Maximal  over  period  vertical  speed’s 
component  for  atmosphere,  heated  by  non-uniform 
wave’s  source  (units  are  dimensionless) 


Figure  12.  Spatial  distribution  of  the  temperature 
for  the  atmosphere,  heated  by  non-uniform  waves’ 
source  (all  units  are  dimensionless) 
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